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In  this  work,  a  transparent  assembly  was  self-designed  and  manufactured  to  perform  ex  situ  experimental 
study  on  the  liquid  water  removal  characteristics  in  PEM  fuel  cell  parallel  flow  channels.  It  was  found 
that  the  dominant  frequency  of  the  pressure  drop  across  the  flow  channels  may  be  utilized  as  an  effective 
diagnostic  tool  for  water  removal.  Peaks  higher  than  1  Hz  in  dominant  frequency  profile  indicated  water 
droplet  removals  at  the  outlet,  whereas  relatively  lower  peaks  (between  0.3  and  0.8  Hz)  corresponded 
to  water  stream  removals.  The  pressure  drop  signal,  although  correlated  with  the  water  removal  at  the 
outlet,  was  readily  influenced  by  the  two  phase  flow  transport  in  channel,  particularly  at  high  air  flow 
rates.  The  real-time  visualization  images  were  presented  to  show  a  typical  water  droplet  removal  process. 
The  findings  suggest  that  dominant  frequency  of  pressure  drop  signal  may  substitute  pressure  drop  as  a 
more  effective  and  reliable  diagnostic  tool  for  water  removal  in  PEM  fuel  cell  flow  channels. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  proton  exchange  membrane  (PEM)  fuel  cell  has  been 
regarded  as  an  ideal  power  source  for  a  variety  of  applications 
due  to  its  significant  advantages,  i.e.,  high  efficiency,  low  emis¬ 
sion,  silence  and  simplicity  [1  ].  The  water  management  in  PEM  fuel 
cell  is  one  of  the  crucial  issues  to  be  fully  understood  and  opti¬ 
mized  before  further  advances  could  be  made.  Sufficient  amount 
of  water  is  necessary  for  maintaining  the  membrane  ion  conduc¬ 
tivity  whereas  excess  water,  or  water  flooding,  may  block  the  pores 
of  electrodes  and  the  flow  channels,  thereby  reducing  the  reactant 
mass  transfer  to  catalytic  sites  [2-5].  The  capacity  to  perform  effec¬ 
tive  water  removal,  therefore,  becomes  one  of  the  most  important 
specifications  for  PEM  fuel  cell  flow  channels  [6].  Since  it  is  highly 
desired  to  know  the  water  behavior  in  both  flow  channels  and 
porous  electrodes,  many  modeling  studies  have  been  conducted 
[7-10].  Several  advanced  considerations  including  two  phase  flow 
in  channels  and  electrodes,  wettability  of  gas  diffusion  layer  (GDL) 
and  transient  processes  in  fuel  cell  have  been  covered  in  these  mod¬ 
eling  studies.  On  the  other  hand,  numerous  experimental  works 
concern  about  the  diagnostic  tools  for  the  water  behavior  in  PEM 
fuel  cell.  Research  in  such  aspect  may  proceed  as  a  more  straightfor¬ 
ward  and  effective  way  to  deal  with  the  water  management,  since 
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the  objective  is  to  sense  or  predict  the  flooding/drying  accurately 
before  applying  corresponding  actions.  Particularly,  transparent 
fuel  cell  was  designed  for  direct  in  situ  visualization  of  water  trans¬ 
port  and  distribution  in  flow  channels  [11-13].  Due  to  the  limited 
view  field  and  water  visualization  effect,  however,  transparent  fuel 
cell  normally  can  apply  straight  parallel  channel  design  only.  In 
addition,  the  water  transport  and  distribution  under  the  land  can¬ 
not  be  revealed  under  transparent  fuel  cell  configuration.  Neutron 
imaging  technique  was  also  utilized  to  obtain  the  water  distribu¬ 
tion  profiles  in  PEM  fuel  cell  [14-16],  with  a  limitation  that  the 
precise  location  of  water  is  difficult  to  determine  due  to  the  over¬ 
lapping  of  images.  Moreover,  magnetic  resonance  imaging  (MRI) 
technique  turned  out  to  be  also  an  effective  tool  to  capture  the 
water  behavior  both  in  membrane  and  flow  field  [  1 7- 1 9 ] .  Its  weak¬ 
ness  is  principally  in  the  requirement  that  the  materials  have  to  be 
nonmagnetic,  making  the  water  in  the  catalyst  layer  and  GDL,  either 
nonwoven  carbon  paper  or  woven  carbon  cloth,  hard  to  visualize. 
A  thorough  review  of  these  techniques  can  be  found  in  Ref.  [20].  In 
addition  to  direct  visualization,  indirect  diagnostic  tool  usually  is  to 
monitor  some  easy-to-measure  engineering  parameters.  Pressure 
drop  across  the  flow  channel  was  proposed  to  be  a  diagnostic  tool 
for  PEM  fuel  cell  flooding  by  Bosco  and  Fronk  in  General  Motors 
[21  ];  similar  findings  concerning  pressure  drop  as  a  diagnostic  tool 
were  reported  by  Barbir  et  al.  [22],  Ma  et  al.  [23]  and  Ito  et  al.  [24]. 
Although  other  parameters,  such  as  cell  resistance  [22]  and  elec¬ 
trode  diffusivity  [25],  were  also  used  as  diagnostic  tool,  pressure 
drop  as  the  mainstream  one  shows  significant  advantages,  namely, 
the  immediate  and  real-time  linkage  between  the  air-water  two 
phase  flow  behavior  and  the  pressure  drop.  Normally,  a  pressure 


1178 


J.  Chen  /  Journal  of  Power  Sources  195  (2010)  11 77- 1181 


drop  increase  is  considered  a  sign  of  water  accumulation  or  water 
film  build-up  in  the  channel;  whereas  a  sudden  decrease  indi¬ 
cates  a  water  removal  at  the  outlet.  However,  pressure  drop  as  a 
diagnostic  tool  for  the  water  behavior,  particularly  water  removal, 
does  have  shortcomings:  first,  high  frequency  oscillation  in  pres¬ 
sure  drop  is  usually  observed,  which  comes  from  the  unstable  two 
phase  flow  transport  in  channel;  such  oscillation  does  not  indi¬ 
cate  any  water  removal/build-up  at  the  outlet  but  make  a  regular 
increase/decrease  pattern  in  pressure  drop  profile  fairly  hard  to 
distinguish;  second,  considering  the  relatively  low  pressure  drop 
(usually  <1000  Pa)  across  fuel  cell  channels,  signal  noises  from  dis¬ 
turbances  and  measurement  errors  may  contribute  considerably  to 
the  results,  reducing  the  accuracy  and  even  changing  the  pattern 
of  the  pressure  drop  profile.  Consequently,  real-time  visualization 
results  are  usually  used  to  assist  the  interpretation  of  the  pressure 
drop  profiles. 

Under  such  circumstance,  the  author  attempted  to  find  a  more 
reliable  and  effective  diagnostic  tool  for  water  removal.  In  perform¬ 
ing  fast  Fourier  transform  (FFT)  upon  certain  amount  of  pressure 
drop  data,  one  can  obtain  the  frequency  spectrum  determined  by 
the  absolute  magnitude  and  phase.  The  normalized  dominant  fre¬ 
quency  is  defined  as  the  weighed  average  of  the  three  components 
with  largest  power  magnitudes.  Compared  with  simply  using  the 
pressure  drop  as  the  diagnostic  tool,  its  dominant  frequency  pro¬ 
vides  more  accurate  and  robust  insights  about  the  water  removal, 
although  the  new  diagnostic  tool  essentially  comes  from  applying 
some  special  mathematical  treatment  on  the  original  pressure  drop 
signal. 

2.  Experimental 

2.1.  Experimental  setup 

In  this  work,  an  assembly  resembling  the  transparent  PEM 
fuel  cell  [11]  was  self-designed  and  precisely  manufactured  using 
CNC  machine.  It  consists  of  an  end  plate  with  two  parallel  micro¬ 
channels  built  on,  a  transparent  plastic  plate  for  visualization  and 
a  window  plate  to  apply  even  compression  to  the  whole  assem¬ 
bly.  Porous  media  (carbon  foam)  was  inserted  into  the  channels 
to  simulate  the  effect  of  GDL  intrusion  in  flow  channel  [26].  Nec¬ 
essary  O-rings  were  placed  on  site  to  avoid  leakage.  The  width  of 
the  channel  is  1.5  mm  and  the  depth  1.2  mm;  consistent  with  typ¬ 
ical  fuel  cell  designs.  The  effective  length  of  the  channel  is  40  mm, 
with  an  entrance  region  to  stabilize  the  air  flow.  Since  the  core 
frictional  pressure  drop  in  the  developing  (entrance)  region  of  the 
channel  is  non-linear  and  unstable  [27],  the  pressure  drop  mea¬ 


surement  was  only  performed  on  the  effective  length  (40  mm)  of 
the  channel.  Pressure  measurement  ports  were  precisely  drilled  at 
two  ends  of  the  channels.  The  experimental  setup  is  schematically 
shown  in  Fig.  1,  with  a  photo  of  the  transparent  assembly  attached. 
The  differential  pressure  across  the  whole  channel  was  measured 
by  a  Setra  230  series  transducer;  data  from  the  transducers  were 
real-time  displayed  and  recorded  using  self-developed  LabVIEW 
programs.  A  CCD  camera  was  attached  to  the  microscope  to  capture 
the  real-time  water  removal  video  at  the  outlet.  Both  de-ionized 
water  and  compressed  air  were  injected  into  the  channels  to  simu¬ 
late  the  two  phase  flow  in  fuel  cell.  Liquid  water  was  injected  into 
each  channel  separately  using  a  syringe  pump,  so  that  the  injection 
rate  can  be  accurately  controlled.  The  water  injection  rate  was  set 
to  0.01  ml  min-1  corresponding  to  liquid  superficial  velocity  {U\)  of 
0.093  m  s-1  in  each  channel,  which  was  an  analogy  to  the  genera¬ 
tion  rate  in  real  fuel  cell  operation  under  a  normal  current  density  of 
0.6  A  cm-2  and  an  effective  MEA  area  of  4  cm2 .  Meanwhile,  air  flow 
rates  were  varied  from  0.01  to  0.81  min-1  using  a  mass  flow  con¬ 
troller,  corresponding  to  stoichiometry  ratios  up  to  10  in  fuel  cell 
operation,  and  gas  superficial  velocities  (Ug)  from  0.09  to  3.7  m  s-1 
in  each  channel. 

2.2.  Mathematical  treatment 

In  this  section,  the  author  would  like  to  present  the  method  to 
obtain  the  dominant  frequency  from  the  pressure  drop  signal.  Sig¬ 
nals  are  converted  from  time  domain  to  frequency  domain  through 
the  Fourier  transform.  It  converts  the  signal  information  to  the 
magnitude  and  phase  component  of  each  frequency.  Customar¬ 
ily  the  Fourier  transform  is  further  processed  to  obtain  the  power 
spectrum,  which  is  the  magnitude  of  each  frequency  component 
squared.  The  spectrum  can  be  then  studied  to  obtain  information 
about  which  frequencies  are  present  in  the  input  signal  and/or 
which  are  significant  component  in  terms  of  the  magnitude  of 
power. 

FFT  is  an  efficient  algorithm  to  compute  the  discrete  Fourier 
transform  (DFT).  Let  x0,. .  .xN_i  be  complex  numbers,  the  DFT  is 
defined  by  [28]: 

JV-l 

Xk  =  YJxne{->2nlN)nk  k  =  0,  . . . ,  k  -  1 

n=0 

Obviously,  this  algorithm  is  complicated  to  work  out  as  it 
involves  many  additions  and  multiplications  of  complex  numbers. 
FFT  is  another  method  for  calculating  the  DFT.  While  it  produces 
the  same  result,  it  is  incredibly  more  efficient,  often  reducing  the 


Fig.  1.  Schematic  of  the  experimental  setup. 
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computation  time  by  hundreds.  The  details  of  FFT  algorithm  are  not 
presented  here  since  it  was  performed  simply  by  relevant  MAT- 
LAB  commands  without  knowing  the  internal  workings.  For  the 
pressure  drop  signal  in  this  work,  it  was  sampled  every  0.16  s  and 
recorded  in  the  data  file.  Therefore  it  can  be  considered  a  discrete 
data  sequence  and  FFT  can  be  performed,  although  the  pressure 
drop  itself  was  analog  signal. 

It  is  necessary  to  introduce  two  self-defined  parameters.  “Win¬ 
dow”  is  the  data  to  be  processed  together  by  FFT.  The  amount  of 
such  data  is  therefore  termed  “window  size”.  One  needs  to  deter¬ 
mine  how  many  data  to  be  processed  at  one  time  (window  size)  so 
that  the  optimum  power  spectrum  that  facilitates  analysis  can  be 
generated.  Details  for  finding  the  appropriate  window  size  are  not 
presented  here  but  can  be  found  in  Ref.  [29].  Another  parameter 
is  defined  to  describe  the  major  frequency  component,  which  is  of 
most  significance  in  the  pressure  drop  signal.  The  power  spectrum 
obtained  after  performing  FFT  on  the  window  will  have  numerous 
frequency  components  that  show  different  power.  The  frequency 
with  the  maximum  power  obviously  is  the  major  component.  How¬ 
ever,  considering  only  the  frequency  with  the  maximum  power 
may  be  inappropriate,  particularly  when  the  second  maximum 
power  is  very  close  to  the  maximum.  Therefore,  the  normalized 
dominant  frequency  is  defined  as  [29]: 

normalized  dominant  frequency 

x  p,  +p2  +  p3  +fl  x  p-1+p2  +  p3  +5  x  Pl  +Pl  +  p3 

where  Pi  is  the  maximum  power  and  f\  the  corresponding  fre¬ 
quency,  P2,  P3  the  second  and  third  maximum  power  and/2,/3  the 
corresponding  frequency,  respectively.  The  normalized  dominant 
frequency  hence  indicated  the  major  frequency  component  of  the 
pressure  drop  signal  by  taking  into  account  the  three  frequencies 
with  largest  power  magnitudes. 

3.  Results  and  discussion 

3.1.  Pressure  drop  profiles 

In  this  section,  the  profiles  of  pressure  drop  across  the  chan¬ 
nels  are  presented  and  analyzed.  Fig.  2  shows  the  real-time  profiles 
when  the  water  was  injected  at  0.1  ml  min-1.  Overall,  higher  air 
supply  rate  resulted  in  higher  pressure  drop.  Apparent  pressure 
drop  oscillations  were  observed  for  air  flow  rate  of  0.05661  min-1 


Fig.  2.  Pressure  drop  profiles  at  water  injection  rate  of  0.01  ml  min-1  and  different 
air  flow  rates. 


Time  (s) 

Fig.  3.  Co-plot  of  dominant  frequency  and  pressure  drop  against  time  at  water 
injection  rate  of  0.01  ml  min-1  and  air  supply  rate  of  0.0566 1  min-1 . 

(the  solid  line).  Every  peak  in  the  oscillation  indicated  one  removal 
of  water  droplet  at  the  outlet,  whereas  the  flat  bottom  of  the  oscil¬ 
lation  corresponded  to  the  injected  water  transport  in  the  channel. 
As  the  water  arrived  at  the  outlet,  it  was  not  expelled  out  imme¬ 
diately  by  the  air  flow;  instead,  it  was  accumulating  to  form  a 
droplet  with  a  diameter  of  ~1.5  mm,  the  droplet  was  not  expelled 
out  until  sufficient  high  pressure  was  gradually  built  up.  The  pres¬ 
sure  drop  sharply  fell  down  as  soon  as  the  droplet  was  removed 
at  the  outlet.  It  should  be  noted  that  although  the  syringe  pump 
ran  continuously  after  the  system  started,  the  actual  water  injec¬ 
tion  via  a  0.7  mm  inner  diameter  hose  was  intermittent,  i.e.,  the 
water  would  not  enter  the  channel  until  sufficient  pressure  build¬ 
up  so  that  the  injection  style  was  droplet  by  droplet.  The  duration 
between  two  peaks  in  pressure  drop  profile  equaled  the  sum  of 
the  waiting  time  between  two  droplet  injections  and  the  time  for 
water  to  transport  through  the  whole  channel.  However,  as  the 
air  supply  rate  was  increased  (the  dashed  and  dotted  lines),  the 
pressure  drop  profiles  became  much  more  coarse  and  irregular, 
making  it  extremely  hard  to  obtain  the  water  removal  informa¬ 
tion  solely  from  the  pressure  drop  profile.  At  increased  air  flow 
rates,  water  could  be  readily  removed  whenever  it  reached  the 
outlet;  the  pressure  build-up/release  cycle  for  water  removal  was 
no  longer  necessary.  As  a  result,  the  water  was  expelled  out  in 
the  form  of  stream,  film  or  tiny  droplet  (stream  being  the  major¬ 
ity),  instead  of  large  droplet  (~1.5mm  diameter)  as  before.  On 
the  other  hand,  the  unstable  pressure  drop  signal  highly  corre¬ 
lated  with  the  two  phase  flow  transport  in  the  channels.  High 
rate  air  flow  was  strongly  disturbed  by  the  water  presenting 
in  channel,  which  explained  the  occurrence  of  many  high  fre¬ 
quency  and  irregular  pressure  drop  oscillations.  For  these  reasons, 
it  may  be  impossible  to  read  the  water  removal  from  the  pres¬ 
sure  drop  oscillation  as  before.  Pressure  drop  as  the  diagnostic 
tool  thus  saw  some  limitation  at  increased  air  flow  rates.  Visual¬ 
ization  images  were  very  necessary  at  such  moment  to  assist  the 
diagnosis. 

3.2.  Real-time  dominant  frequency  plots 

The  normalized  dominant  frequency  and  the  original  pressure 
drop  were  co-plotted  against  time  in  Fig.  3  when  the  water  injection 
rate  was  0.01  ml  min-1  and  air  supply  rate  0.0566 1  min-1 .  To  obtain 
the  frequency  profile,  at  any  given  time  except  the  first  second,  the 
window  was  applied  on  certain  amount  (window  size,  in  this  case 
36)  of  data  right  before  the  datum  at  that  moment.  Therefore,  the 
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advance  of  frequency  change  to  diagnose  the  water  behavior  could 
be  guaranteed.  Note  that  at  the  very  beginning  (the  first  second) 
of  data  recording,  the  window  was  applied  on  the  data  right  after 
since  no  historical  data  available  yet. 

The  water  removal  pattern  in  this  case  may  be  named  “inter¬ 
mittent  droplets  mode”.  As  clearly  illustrated  in  Fig.  3,  every  time 
before  the  water  removal,  the  dominant  frequency  experienced 
a  peak  together  with  the  pressure  drop.  Also,  the  real-time  visu¬ 
alization  indicated  that  the  magnitude  of  the  peak  in  frequency 
correlated  with  the  size  of  the  water  droplet  that  was  removed 
at  the  outlet.  For  example,  the  ~1.8Hz  peak  at  about  80s  corre¬ 
sponded  to  the  removal  of  a  ~1.5mm  diameter  droplet;  whereas 
lower  peak  at  about  25  s  to  smaller  water  droplet  with  diameter 
of  ~  1.0  mm.  It  is  worthy  to  note  the  last  two  peaks  in  dominant 
frequency  locating  within  the  peaks  of  pressure  drop,  which  are 
much  wider  than  previous  ones.  The  visualization  indicated  that 
during  the  wider  peak  in  pressure  drop,  the  water  was  accumulat¬ 
ing  and  then  forming  droplet  but  differently,  it  was  not  expelled  out 
as  the  droplet  was  formed;  the  droplet  grew  larger  and  eventually 
broken  down  into  a  series  of  small  droplets  (diameter  ~0.3  mm), 
which  then  flowed  out  very  quickly.  The  unique  peak  in  domi¬ 
nant  frequency  during  such  process  perfectly  indicated  the  one 
removal  of  a  series  of  small  droplets  and  predicted  the  moment 
of  removal  (right  after  each  peak).  Flowever,  judging  from  the 
pressure  drop  profile,  one  may  get  confused  at  the  wider  peaks 
and  incorrectly  speculate  that  the  water  was  removing  during  the 
whole  peak. 

A  typical  water  removal  process  is  presented  in  Fig.  4.  The  still 
images  can  only  show  part  of  the  whole  process  with  the  critical 
scenes.  Specific  time  for  taking  each  scene  is  indicated  so  that  a 
direct  comparison  can  be  made  with  profiles  in  Fig.  3.  From  Fig.  4(a) 
to  (c),  the  pressure  was  gradually  built  up  with  water  accumulation 
and  droplet  growth,  during  which  the  pressure  drop  profile  expe¬ 
rienced  apparent  uptrend.  The  edge  of  the  formed  droplet  could  be 
clearly  seen  in  these  figures.  In  Fig.  4(a),  the  water  from  the  channels 
accumulated  at  the  outlet;  then  a  droplet  was  formed  and  began  to 
vibrate  (the  unstable  edges  could  be  clearly  seen  in  (b)  and  (c)).  In 
Fig.  4(d)  and  (e),  the  water  was  promptly  flowing  out  in  droplet, 
corresponding  to  the  abrupt  pressure  drop  decrease  around  180  s 
in  Fig.  3. 

Fig.  5  shows  the  dominant  frequency  and  pressure  drop  co¬ 
plot  when  the  air  supply  rate  was  increased  to  0.23661  min-1.  In 
this  case,  the  water  removal  at  the  outlet  was  not  so  regular  like 
previous  case;  the  droplet  formation  and  growing  process  were 
not  frequently  observed.  Instead,  the  water  was  flowing  out  in 
stream  for  most  of  time,  corresponding  to  many  relatively  lower 
dominant  frequency  peaks  (0.4-0.6  FIz)  in  Fig.  5.  The  dominant  fre¬ 
quency  peak  higher  than  1 .2  Hz  indicated  there  was  still  one  water 
removal  in  droplet  (~1 .0  mm  diameter)  at  around  1 20  s,  consistent 
with  the  visualization  result.  The  increased  air  supply  rate  obvi¬ 
ously  made  the  water  removal  more  chaos,  as  the  increased  air 
flow  could  readily  provide  the  pressure  needed  to  push  out  the 
water  at  the  outlet.  Consequently,  instead  of  the  intermittent  mode 
in  previous  case,  the  water  removal  could  be  considered  almost 
continuous,  not  completely  because  there  were  still  low  frequency 
regimes  (<0.3  Hz)  in  Fig.  5.  Examining  the  pressure  drop  profile, 
one  cannot  obtain  the  water  removal  information,  since  the  pres¬ 
sure  drop  was  influenced  by  the  water  transport  in  channel  also; 
particularly  at  increased  air  flow  rate,  the  instability  effects  were 
amplified.  Fig.  6  presents  the  visualization  image  for  water  trans¬ 
port  in  channel  at  air  supply  rate  of  0.2366 1  min-1 ,  which  may  help 
the  understanding  of  the  unstable  pressure  drop.  It  can  be  seen 
that  water  transport  occurred  both  in  forming  a  film  at  the  channel 
sides  and  in  passing  through  the  pores  of  porous  media  (a  represen¬ 
tative  region  shown  in  the  ellipse).  Such  water  transport  in  channel, 
no  matter  in  either  way,  would  disturb  the  air  flow  and  consider- 


Fig.  4.  Visualization  images  of  a  typical  water  removal  process  occurred  at  water 
injection  rate  of  0.01  ml  min-1  and  air  supply  rate  of  0.0566 1  min-1 :  (a)  water 
accumulation  at  t=  173.2  s;  (b)  droplet  formation  at  t=  175.3  s;  (c)  pressure  further 
build-up  and  droplet  being  pushed  at  t  =  1 76.8  s;  (d)  the  beginning  of  droplet  removal 
at  t=  178.3  s;  (e)  water  removal  and  pressure  release  at  t=  179.5  s. 

ably  contribute  to  the  instability  of  pressure  drop  profile.  Therefore, 
although  one  could  find  many  “peaks”  in  the  pressure  drop  pro¬ 
file  in  Fig.  5,  most  of  them  were  not  from  the  water  removal  at  the 
outlet  but  the  unstable  two  phase  flow,  particularly  the  water  trans¬ 
port  in  the  channel.  However  for  the  frequency  profile,  the  lower 
peaks  (0.4-0.6  Hz)  in  dominant  frequency  indicated  water  removal 
at  the  outlet  in  form  of  stream  and  higher  peak  (>1.2  Hz)  in  form 
of  droplet.  The  dominant  frequency  profile  clearly  shows  advan¬ 
tages  over  pressure  drop  profile  in  diagnosis  of  water  removal  at 
the  outlet. 
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Fig.  5.  Co-plot  of  dominant  frequency  and  pressure  drop  against  time  at  water 
injection  rate  of  0.01  ml  min-1  and  air  supply  rate  of  0.2366 1  min-1 . 


Fig.  6.  Visualization  image  of  the  water  transport  in  channel  at  water  injection  rate 
of  0.01  ml  min-1  and  air  supply  rate  of  0.2366 1  min-1 . 

Fig.  7  shows  the  co-plot  obtained  at  further  increased  air  sup¬ 
ply  rate  0.32661  min-1.  The  pressure  drop  profile  contains  more 
intense  oscillations  as  the  water  transport  in  channel  became  more 
unstable.  However,  applying  the  above  criteria  on  the  dominant  fre¬ 
quency  profile,  one  can  still  find  that  the  water  droplet  removals 
took  place  at  approximately  15,  45,  55,  85  and  185  s,  and  several 
water  stream  removals  at  other  moments.  Again,  the  dominant  fre¬ 
quency  as  a  diagnostic  tool  shows  better  reliability  and  resistance 
to  noise  than  pressure  drop,  particularly  at  increased  air  flow  rates. 


Fig.  7.  Co-plot  of  dominant  frequency  and  pressure  drop  against  time  at  water 
injection  rate  of  0.01  ml  min-1  and  air  supply  rate  of  0.3266 1  min-1 . 


4.  Conclusions 

In  this  study,  the  relations  between  the  water  removal  in  the 
PEM  fuel  cell  parallel  channels  and  the  dominant  frequency  of 
pressure  drop  signal  were  experimentally  investigated  using  a  self- 
designed  and  manufactured  transparent  assembly.  It  was  found 
that  dominant  frequency  showed  better  performance  in  water 
removal  diagnosis  compared  with  pressure  drop,  which  presently 
is  the  principal  diagnostic  tool  for  the  water  behavior.  The  real¬ 
time  visualization  results  indicated  that  peaks  higher  than  1  Hz  in 
dominant  frequency  corresponded  to  the  water  droplet  removals 
(after  its  formation  and  growth)  and  peaks  between  0.3  and  0.8  Hz 
corresponded  to  the  water  stream  removals.  Additionally,  higher 
magnitude  in  dominant  frequency  peak  corresponded  to  water 
droplet  with  larger  size.  The  pressure  drop  across  the  channels 
may  not  diagnose  the  water  removal  accurately  because  it  was 
readily  influenced  by  the  two  phase  flow  transport  in  channel, 
although  also  correlated  with  the  water  removal  at  the  outlet.  Par¬ 
ticularly  at  high  air  low  rates,  dominant  frequency  showed  much 
better  reliability  and  higher  accuracy  than  pressure  drop  for  water 
removal  diagnosis.  It  may  be  preferred  to  use  dominant  frequency 
as  an  effective  diagnostic  tool  for  water  removal  in  PEM  fuel  cells. 
Future  work  will  mainly  concern  the  establishment  of  a  math¬ 
ematical  model  correlating  the  dominant  frequency  and  water 
removal  characteristics  such  as  droplet  diameter  and  cycle  dura¬ 
tion. 
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